Atom-molecule coherence in a one-dimensional system 
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We study a model of one-dimensional fermionic atoms with a narrow Feshbach resonance that 
allows them to bind in pairs to form bosonic molecules. We show that at low energy, a coherence 
develops between the molecule and fermion Luttinger liquids. At the same time, a gap opens in the 
spin excitation spectrum. The coherence implies that the order parameters for the molecular Bose- 
Einstein Condensation and the atomic BCS pairing become identical. Moreover, both bosonic and 
fermionic charge density wave correlations decay exponentially, in contrast with a usual Luttinger 
liquid. We exhibit a Luther-Emery point where the systems can be described in terms of nonin- 
teracting pseudofermions. At this point, we provide closed form expressions for the density-density 
response functions. 
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Recent experiments on fermionic 6 Li or 40 K atoms in 
optical traps have led to the realization of paired super- 
fluidity. Pairs of atoms were found to bind into bosonic 
molecules that displayed a Bose-Einstein condensation^ 
0, E3- H as the magnetic field was varied across a Fesh- 
bach resonance. The Bose condensation of the molecules 
is expected to trigger superfluidity of the fermions. A 
crossover is thus expected from a BCS superfluid of 
paired atoms to a Bose Einstein condensate (BEC) of 
molecules as the molecules become more tightly bound. 
It is important to stress that since the BEC and the 
BCS states have the same broken symmetry, there is no 
fundamental distinction between them, and these two 
extreme limits are connected by a smooth crossover. 
This crossover is naturally described by the fermion- 
boson model which has extensively been considered in the 
context of bipolaronic and high-Tc superconductivity 5] 
and has recently known a regain of theoretical interest 

|M 1^1 - 

Recently, the BCS-BEC crossover has been investi- 
gated in the one-dimensional case0, Q, 0|. Despite 
the absence of broken symmetries in one dimension, a rich 
phenomenology is known to emerge such as collectiviza- 
tion of single-particle degrees of freedom, spin-charge sep- 
aration and quasi-long range order^^. In particular, 
spin-1/2 fermions with attractive interactions give rise 
to a state with gapless charge degrees of freedom and 
gapped spin degrees of freedom, with quasi-long range 
superconducting and char ge d ensity wave order, known 
as the Luther- Emery liauidjlflj. In [ 1 it 1 1 2l | . an integrable 
modeljlj] of spin-1/2 fermions with attractive interac- 
tions was considered. It was shown that the Luther- 
Emery liquid obtained at small attraction was cross- 
ing over to a Luttinger liquid of tightly bound bosonic 
molecules. In 10], the boson- fermion model was used, 
but the case of a broad resonance where only the bosons 
or only the fermions are present was considered, yielding 



results similar to those of the integrable model. An in- 
teresting theoretical question is to understand what hap- 
pens in the case of a narrow Feshbach resonance when the 
atoms and the molecules can coexist. In the present Let- 
ter, we investigate such a case. We show that a strongly 
correlated state exists in this case, in which the order pa- 
rameters of the Bose condensation and superfluidity de- 
cay with the same critical exponent, and density fluctua- 
tions near the Fermi wavevector are strongly suppressed. 
Also, a spectral gap is present in the atomic spectral 
function, as in the case of a Luther-Emery liquid. From 
the experimental point of view, the interest in working in 
one dimension is the possibility of using the confinement 
induced resonance (CIR) to form the molecules, as 
recently demonstrated 1 1 fij . Another way of forming the 
molecules is by photoassociation techniques |ljj. Since 
it is possible to realize quasi-ld systems of bosons and 
fermions in optical traps [lq. some of our predictions 
could be testable in experiments. 

We consider a ID system of fermionic atoms that can 
bind reversibly into bosonic molecules. In the continuum 
case, the Hamiltonian of the system reads: 
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dxdx'V vv >(x - x')p v (x)p v >{x') (1) 



where ip a annihilates a fermion (atom) of spin a, tpb an- 
nihilates a boson (molecule), p^ — iplipb Pf — ^2 a 
and we have set % = 1. The parameters Vff, Vbb and 
Vbf measure (respectively) the fermion- fermion, boson- 
boson, and fermion-boson repulsion. The parameter v 
is the detuning, and the term A describes the binding of 
pair of atoms into a molecule and the reverse process. 
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When A 7^ 0, only the total atom number AT = 2Nb + Nf 
is conserved. The fermion-boson model can be ana- 
lyzed by bosonization techniques in the limit of a narrow 
resonance (i.e. small A) and provided that neither the 
density of atoms nor the density of molecules vanishes. 

First, we recall the bosonized description of the sys- 
tem when boson-fcrmion interactions are turned off. For 
A = 0, Vbf = 0, both the number of free atoms Nf and 
the number of molecules Nb are conserved and the Hamil- 
tonian equivalent to Q is given bvfliij]: 

2<?i± 



H 



Hb + H n 
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dx cos v^ct (2) 



dx 
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where [<f> v (x), n„/ (x 1 )] = iS(x — x')5 v y , (y, v' — b,a,p). 
The parameters u p , u a , Ub are the velocity of respectively 
atomic density, atomic spin, and molecular density ex- 
citations. K p and Kb are the Luttinger exponents 0|. 
They decrease as (respectively) the atom-atom and 
molecule-molecule repulsion increase. For weak repul- 
sion, it is possible to express K p as a function of the 
scattering length a Sl trapping frequency lo± and the 
Fermi velocity of non-interacting atoms as K p = (1 + 
2a s uj±/irvp)~ 1 / 2 . In the case of bosons, the Luttinger 
exponent Kb must be extracted from the Lieb-Linigcr 
equations [Lj,ll9|. For weak interaction, Kb — * +oo, while 
in the Tonks- Girardeau limit, Kb = 1. Turning to the 
spin interaction, it is known 13J that under the renoma- 
lization group (RG) gi± and K a flow to the fixed point 
K* = 1 , glj_ = provided that the repulsive interactions 
respect the SU(2) spin symmetry. We will thus replace 
Ka,gi± by their fixed point value in the following. The 
power of bosonization in treating interacting systems in 
one-dimension comes from the possibility of expressing 
the fermion and boson annihilation operators in terms of 
the fields in @ as: 



e j S Wp-r<t>p+a(0<r-r<t>< r )Kx) 



V2t 



(4) 



[1+Acob(206-2* b x)], (5) 



where the dual fields^] v {x) — ir j x Ii v (x')dx' (y = 
b,p,cr) have been introduced, kp = nNf/L, ks — 
irNb/L, r = ±, and a is a cutoff. The atom field is 
expressed in terms of the right and left moving fields of 



n — ik 



FX ' l i ) -,i^. Similarly, 



S^-i' 



Eq. as: ij) a {x) = e ikFX ip +:C 
the atom and molecule density are given by 

s/2 cos(2k F x - V2<j) p ) 

Pf{x) = O x (j)p-\ ^cos 

7r ira 

1 co$,(2k B x - 2<f) b ) 

p b {x) = d x <j>b H . 

7r ira 

Let us now turn on A, Vbf,v ^ 0, assuming that that 
they are small compared with the kinetic energy of 



V2<M6) 
(7) 



the atoms and the molecules. This corresponds to the 
limit of a narrow Feshbach resonance [13. For k B ^ 
kp, the boson- fermion repulsion reduces to a term ~ 
/ d x (j) b d x (f> P ^^- Using Eqs. Q-© The bosonized 
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form of the A is: 
tibf - 



dx cos(6b — V28 p ) cos \f24> a (8) 



Finally, the detuning term reads: — ^ J dxd x <f>b, which 
shows that is can be eliminated by a shift of cf>b — > 
<t> b - 2^|ps without affecting Eq. Thus, the de- 

tuning is effective only when it induces band-filling tran- 
sitions at which either the density of atoms [y < 0) or 
of molecules {y > 0) vanishes. At these transitions, the 
bosonization descri ptio n breaks down, with divergence of 
the compressibility. [21| The large detuning limit has been 
analyzed in [Toj ]. where it was shown that existence of 
virtual atom or molecule states only leads to a renormal- 
ization of the Luttinger exponent Kb or K p respectively. 
In the following, we take v = 0. The RG equation for A 
reads: 
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showing that for ^ jj^- < 3/2, this interaction is 

relevant. Since for hard core bosons Kb = 1 and for 
free bosons Kb = oo while for free fermions K p — 1, this 
inequality is satisfied except for very strongly repulsive 
interactions Vff and Vbb- In this strong repulsion limit, 
there is no coherence between atoms and molecules and 
the system is described by the theory of |20|. For less 
repulsive interactions, the relevance of A drives the sys- 
tem to a new fixed point. To understand the nature of 
this fixed point, it is convenient to perform a canonical 
transformation 13] 0_ = ^0 b -^6 p , 9+ = ^d b + ^0 p 
and the same transformation for the cj> v . Then, Hb + H„ 
becomes H± with: 



H+ = 



dx \ -> 

9^ 2-^, 



L w K v {^f + -^{d x ^f 



f dx 

+ J — [5i(7rn + )(7rn_)+ 52 9 x + 9 x 0_],(lO) 

where u±, K±, <?i.2 can be obtained in terms of the pa- 
rameters of the original Hamiltonian and VJ,/. In the 
following, we will assume u± — u a — u. When we ex- 
press (JBJl in terms of we see that when A is rele- 
vant, only the fields CT and are locked, with a gap 

A_ = A CT ~ «/«(Aa/«) 4/(6 " 2 ^ 1 "^ 1) , while 0+ re- 
mains gapless. The field </>+ describes the total density 
excitations of the system, as can be seen from the rela- 
tion \/6[0+(oo) — (f> + (— oo)] = TV. It is described by the 
low-energy Hamiltonian: 



dx 
2^ 



u* + K* + (nU + ) 



K* 



■(5^+) 2 



(11) 
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where u*,, K5_ are renormalized values of u + , K + , result- 
ing from the residual interactions with the gapped mode 
caused by the terms <?i,<?2- The gapful excitations 
are formed by kinks of the fields 6>_ and such that 
6- (oo) — d- (—00) = ±n/^/3 and <p a+ (00) — <p a _ (—00) = 
±n/>/2. As a result, they carry a spin ±1/2 and create 
a phase difference of ^fir between bosons and fermions. 
Thus, they can be identified with half-vortices carrying a 
spin 1/2. Having understood the spectrum of the system, 
let us turn to its ground state correlations. The locking 
of the fields 6L and CT yields from Eq. (J5J the follow- 
ing low energy expression for the BEC order parameter 
for the molecules: b{%) ~ e 4 V / ¥ e +, while the order pa- 
rameter for s-wave BCS superfluidity of the fermions be- 
comes: VtV'J. ~ e lv/Ie 'cosV2^ ~ eVP+. Thus, these 
two order parameters become identical in the low en- 
ergy limit as in the case of higher dimensionality. 0, . 

The boson correlator behaves as: (^^(x, t)4 , b(0 j 0)) = 

1 

((x 2 + u 2 T 2 )/a 2 ) 6K + , yielding a molecule momentum 
distribution Tis(fc) ~ l^l 1 '' 3 -^) -1 . This momentum dis- 
tribution could be measured in a condensate expansion 
experiment |23l l24j|. A more striking consequence of the 
locking 9_ is that both the ±2kg and the ±2kp harmonic 
in (respectively) pb(x) and Pf{x) acquire exponentially 
decaying correlations, with a correlation length ~ u/A_. 
The origin of this exponential decay is that both pb and 
Pf depend on the field (j)- dual to Such behavior 

is in contrast with the behavior of a Luttinger liquid of 
molecules, or a Luther Emery liquid of fermions with at- 
tractive interaction, in which these correlation functions 



would decay as power law. 

A more detailed picture of the gapful spectrum and 
the correlation functions can be obtained at a particu- 
lar solvable point of the parameter space, the so called 
Luther-Emery point |25|. For K- = 3/2, one can intro- 
duce the pseudofermion fields: 



*[(vl^-- r vi e -)+ CT ( e -- r ^)] 



V2t 



(12) 



to rewrite Eq. © and H- + H a as a free pseudofermions 
Hamiltonian: 



H = J dx 



~iu ^ r^l i(T d x ^ r ,a 
r=± 



A 



(13) 



At that point, the kinks are pseudofermions with dis- 
persion e(k) = zb-y 7 '(uk) 2 + A 2 , where A = ^§^< In 
terms of the pseudofermions the fermion density reads: 



H. c. which yields an expression of density-density cor- 
relations correlation in terms of modified Bessel func- 
tions. The boson density is expressed in terms of 
order and disorder operators of four 2D non critical 
Ising models^ as: p b {x) ~ e i{y/m<t>+-^ B x)^ i + 

«o'i,TA t 2,T)(/ i i,l cr 2,i +* cr i.iA t 2.i) + H. c. where <j{p) are the 
order (disorder) parameters of the Ising model. The cor- 
relation functions of the non-critical Ising model are ob- 
tained from [23. The Fourier transform of the Matsubara 
correlation functions then reads: 
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and the response functions are then obtained by the 
substitution ilu — > u> + iO. Since the generalized hy- 
pergeometric functions p +iF p (. ..;...; z) are analytic for 
\z\ < 1 [23, the imaginary part of the response func- 
tions vanishes if oj < 2A. For li > 2A, the behavior 
of the imaginary part is given by the expression J29| of 
the imaginary part of 3.F2 in terms of Appell's hyperge- 



ometric function[2|| F3, yielding power law singularities 
at u> = 2 A with an exponent -K+/3 — 1/2 for the atoms 
and — 1/2 for the molecules. The Luther-Emery 

limit also yieds the q ~ components of the density re- 
sponse. Noticing that ps — d x <j)+/V3 + \/2d x <f>-/V3, we 
find that 3Xp P ((j,w) is the sum of a term cx 6(lj ± u + q) 
coming from the gapless phase mode, and a term cx 
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(w - 2A - (wg) 2 /(4A 2 ))- 1 /2 e ( ti , _ 2A - (uq) 2 /(4A 2 )) 
coming from the gapped mode. The same result holds 
for 3xJ p (q,[j). Most interestingly, the cross correla- 
tions of the fermion and the boson density are also non- 
vanishing for q — > and behave similarly to x^ p (g,w). 
The imaginary parts of correlation functions Eq. I|14f> - 
(|15|1 can be measured by Bragg spectroscopy with large 
momentum transfer |30j| and the q = component can be 
measured b y B r agg spectroscopy with small momentum 
transfer [sil l33j . The atom Green's function is ob- 
tained using form factor expansion techniques |34ll35| as: 



K + < 1/3 i.e. very strong repulsion. In the Mott in- 
sulating state, superfiuid and BEC correlations become 
short ranged as the density density correlations. 

Note added: The model considered in the present pa- 
per has also been studied independently by D. E. Sheehy 
and L. Radzihovsky 3SJ . The CSG phase discussed by 
these authors is identical to the one discussed in the 
present paper, and the condition for stability of the CSG 
phase derived by them is also compatible with the one 
derived in this paper. 



G(x,T)~e«<> - 
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where p = \J x 1 + v 2 r 2 .The corresponding spectral func- 
tion is obtained by Fourier transformation]^ .and it 
vanishes below the gap A, as in a superfluid|36|. It 
would be interesting to observe this gap in a conden- 
sate expansion experiment |U |2^. We have seen that 
with a narrow resonance, the mutual coherence of the 
atoms and the molecules reinforces the superfluidity of 
the system. An important question to ask is how such 
a coherence can be lost. A first way of losing the co- 
herence is by applying a temperature, creating a density 
of half-vortices ~ e~ A "/ T , which destroys phase coher- 
ence between the atoms an molecules on a lengthscale 
£(T) - e A "/ T . The second way of losing the coherence 
between bosons and fermions is by applying a magnetic 
field strong enough to cancel the gap for the creation of 
half- vortices excitations, causing a commensurate incom- 
mensurate transition. |37j As a result of this transition, 
power law singularities in the density-density correlations 
reappear, and the behavior of the system is again de- 
scribed by the models of Refs. [2(J. ^From the experimen- 
tal point of view, a narrow resonance could be obtain by 
working with 6 Li atoms0 trapped in a two dimensional 
optical lattice. The relevant parameters 0? ] for 6 Li are 
the mass m( 6 Li)=9.96 x 10~ 27 kg, the width of the reso- 
nance AB = 0.23G = 2.3 x KT 5 T, the atom-atom scat- 
tering length db g =80 do = 4.23 x 10~ 9 m, the difference 
in magnetic moment between the atom and the molecule, 
Afi ~ 2/i B =2 x 927.400949 x 10~ 26 = 1.8 x 10~ 23 
J.T -1 , and the trapping frequencies lu±~2tt x 69kHz, 
uj z = uj±/ '270. With these parameters we find that the 
Fermi velocity is of the order 10 -2 m/s and the value of 
the Luttinger parameter is K p ~ 0.995. Finally we would 
like to comment on the fact we have only considered the 
case of the continuum system JJJ . All the considerations 
above apply equally to a lattice model at an incommensu- 
rate filling, albeit the effective mass of the atoms and the 
molecules can be strongly enhanced by the periodic po- 
tential. For commensurate filling in the lattice system, an 
umklapp term, H urn k. — / dx cos \/2A(j) + is present 

in the Hamiltonian and can create a Mott gap0|. The 
RG treatment shows that the Mott gap exists only for 
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